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TECHNICAL PAPER 


FORBIDDEN TANGENTIAL ORBIT TRANSFERS BETWEEN 
INTERSECTING KEPLERIAN ORBITS 

INTRODUCTION 


Although the subject of impulsive transfers between Keplerian orbits has been treated in many 
papers (e.g. [ 1 ,2]), there has been no demonstration that certain transfers are not allowed. In this paper, it is 
shown that, for planar intersecting orbits, the tangential impulse transfer cannot always be utilized. The 
treatment is formulated as purely geometric in nature and no consideration is given to optimality. 

For the two-impulse transfer between two coplanar elliptical orbits, there exists a triple infinity of 
elliptic orbits connecting the given orbits. In the case of realistic space transfer maneuvers, the optimum 
impulses are contained in a small useful angle between the tangential direction and the local horizontal 
direction as seen on the transfer orbit [3], As a practical application, we may search for the optimum 
cotangential transfer as a sub-optimum solution. 

The tangency conditions reduce the triple infinity of solutions to a single infinity of solutions so we 
are at liberty to parameterize these solutions by one convenient variable; the true anomaly of the initial orbit 
is chosen. It is then demonstrated that the descriptors of the transfer orbit can either be determined directly 
in terms of this parameter or eliminated from the discussion. Previous studies relied on a series of sequential 
calculations, and this obscured the fact that singularities can sometimes occur in the semi-major axis and 
eccentricities of the transfer orbit. 

The arguments which exclude certain transfers lor intersecting orbits are shown to be inapplicable 
for the case of non-intersecting orbits; singularities occur only for transfers between intersecting orbits. 

I thank the (anonymous) reviewer who pointed out an interesting geometric argument based on the 
work of Professor Busemann [4] for Reference 5. The comments that were returned to me used a Busemann 
space construction to show that if the initial and final orbits do not intersect then the transfer solution is 
always available; in the case of intersecting orbits, forbidden transfer regions can occur. This entirely 
agrees with the results presented here. 

The exact conditions for the occurrence of forbidden transfers are presented in this paper. Their 
existence provides a new result in a classical subject. 


PROBLEM FORMULATION 


In the context of the (planar) Keplerian two body problem, define a cartesian coordinate system 
with the attracting primary at the origin. The initial orbit is oriented on this coordinate system in such a way 
that its perigee occurs along the x-axis. The initial orbit will be taken to have an eccentricity ei and a semi- 
major axis a i . The final orbit has corresponding parameters e 3 and a 3 and is oriented in such a way that its 
semi-major axis is rotated clockwise through an angle ft from the negative x-axis. See Figure 1 . 



An orbit transfer from orbit No. 1 to orbit No. 3 can occur in many ways but only a subset of general 
transfers will be considered here, namely tangential transfers. The geometry of a transfer orbit can now be 
taken as equivalent to constructing a confocal conic section from some point, <J>, , on the initial orbit to a 
point <j >2 on the target orbit. 

Four constraints exist, i.e., the radius of the transfer orbit must match at the initial and final points 
and the velocities must be colinear at the same points. The unknowns of the problem include the semi-major 
axis and eccentricity of the transfer orbit, the rotation (phase) angle of the transfer orbit, the angle 4> i (the 
angular leave point) and the angle 4> 2 (the angular arrival point). Since there are five unknowns and only 
four constraints, it is apparent that one free parameter exists and it is reasonable to parameterize with 
respect to the leave angle, <J>|. Figure 2 shows a typical transfer orbit. 
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Figure 2. 

To proceed with a quantitative formulation, the three orbits can be characterized by 

r = a, (1 -cr) ( 1 ) 

1 1 -t- e i cos c{> 


a 2 (1 - e 2 2 ) 

1 + e 2 cos(4> + p) 


( 2 ) 


— a 3 ( * ~ e 3~) (3) 

1 + e 3 cos(tJ> + i)) 

where T| , r 2 , r 3 are the radii of the three orbits, a, , a 2 , a 3 and e, , e 2 , e 3 are the semi-major axes and eccentrici- 
ties of those orbits. The angle p is the rotation angle of the apogee of the transfer orbit, is the rotation 
angle of the apogee of the final orbit, and (j> is the true anomaly of the initial orbit. 

From these 


dr, _ r, 2 e, sin <j» 

d4> a | (I — e i *") 


(4) 
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dr 2 

fi(J> 


r 2 2 e 2 sin(4> + p) 
a 2 ( 1 - e 2 2 ) 


( 5 ) 


dr 3 r 3 ~ e 3 sin(<t> + v)) 
dct> a 3 (l-er) 


Equating the radii and requiring tangency at a generic point <J>, gives the first two constraint equations 
a, (1 -e, 2 ) a, ( 1 -e-> 2 ) 

: (7) 

l+e|Cos<p| 1 +e-> cos(4>| + p) 


r i 2 ei sin (J>, _ r 2 2 e 2 sin(4> j + p) 
a,(I-e, 2 ) a 2 ( 1 - e 2 2 ) 


( 8 ) 


and the second two constaint equations come from a similar generic point cj> 2 , i.e.. 


a 2 ( I e 2 ~) a 2 ( 1 e 3 ~) 

1 + e 2 cos(<J) 2 + p) I + e 3 cos(4» 2 + 0) 


(9) 


r 2 2 e 2 sin(4> 2 + p) 
a 2 ( I — e 2 2 ) 


r 3 2 e 3 sin(<{> 2 + 0) 
a 3 ( 1 - e 3 2 ) 


( 10 ) 


Since r ( (4> , ) = r 2 ((J)|) and r 2 (cj) 2 ) = r 3 (4> 2 ) and abbreviating 


Q, 


I 

aid -er) 


(ID 


equations (7) through (10) can be written as 


Qi( I + ei cos (J)| ) = Q 2 [ 1 + e 2 cos(<j>i + p)J 


( 12 ) 
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Q, C| sin 4>| = Q 2 e 2 sin(<|>, + p) 


(13) 


Q 2 [ 1 + e 2 cos(c(> 2 + p)] = Q 3 [ 1 + e 3 cos(<j) 2 + -9)] 


(14) 


Q 2 e 2 sin(4» 2 + p) = Q 3 e 3 sin(<t> 2 + ft) 


(15) 


It should be noted that equations ( 1 2) through ( 1 5) constitute four equations in five unknowns, i.e. , 
Q 2 , e 2 , p, c{> 2 and the parameter <j>|. 

The first solution to those equations which were obtained was the most general case, i.e., e, 0 
and e 2 t 0 Subsequently, it was found that the singular cases where either e, ore 3 (ore, ande 3 )vamshe 
should be treated separately. Since the derivation for singular cases is much easier than the general case the 
order in which the solutions were obtained will be reversed, i.e. , the presentation wi sta W1 e ' J* 

0 (Hohmann transfer) , then the two cases where either e , = 0 or e 3 = 0, and finally the genera so u i 
will be seen, later, that the general solution does not readily reduce to the special cases. 

First Singular Case: Hohmann Transfer 

Setting e, = e, = 0, and. obviously. 9 = 0, in equations (12) through (15) gives 


Ql = Q 2 [l + e 2 cos(<J>! + p)] 


(16) 


O = Q 2 e 2 sin(cf>i + p) 


(17) 


Q 2 [l +e 2 cos(4> 2 + p)] — Q 3 


(18) 


(19) 

Q 2 e 2 sin(4> 2 + p) — 0 . 


Excluding Q 2 = 0, e 2 = 0, and Q, = Q 3 immediately shows that, from (17) and (19) 


COS(4> 1 + p) — — 1 


COS((t> 2 + p) — ± 1 


5 



which leaves the two equations 


Qi = Q2 (I ±e 2 ) 

Qi = Q: (1 ±e 2 ) 

Since Qi + Q 3 the immediate conclusion is that there are two cases. Either 

Qi = Q2 (1 + e 2 ) 

and 

Q? = Q2 ( I — e 2 ) 

or 

Qi — Q2 (1 ~e 2 ) 

and 

Q3 = Q 2 ( 1 + e 2 ) 

From the first case comes the solution 

Q. -Q.3 

62 Q.3 + Q. 


Q 2 = 


Qi + Q3 

2 


(20) 


( 21 ) 


Since e 2 > 0, this solution applies whenever Q, > Q 3 . From the second case. 


= Q3-Q1 
62 Qj + Qi 


(22) 
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_ Qi + Qj 

Qi - 2 


( 23 ) 


which should be used when Q 3 > Q|. 

4>i is arbitrary, and for the transfer to be meaningful, p = ±tt. Then, 4> 2 = ± tt - 

Second Singular Case: Initial Eccentricity = 0, Final Eccentricity + 0 

Since e, = 0 the only thing which vitiates the spatial isotropy is the eccentricity of the target orbit. 
The rotation angle of the apogee of the target orbit can then be chosen to be equal to 0 (i.e. , ft = 0) and 
equations (12) through (15) can be written as 


Qi = Q 2 [l + e 2 cos(ch + p)] < 24) 

0 = Q 2 e 2 sin(4>, + p) (25) 

Q 2 [l + e 2 cos(4> 2 + p)] = Q 3 [l + e 3 cos (j> 2 ] (26) 

Q 2 e 2 sin(cf) 2 + p) = Q 3 e 3 sin 4> 2 (27) 

From equation (25), either 

p=-4>, (28) 

or 

p = TT — 4> | (29) 


so two choices exist for p. Using a “parallel” derivation to reduce space produces the following. 


7 


If p = — 4>i then 


If p = tt — 4>i then 


Qi = Q:(l + e 2 ) 

(30) 

Qi ~ Q 2 ( 1 — c 2 ) 

(30) 

Q:1 1 + cos( 4 > 2 - <J>i)] = Qd 1 + e 3 cos<)) 2 ] 

(31) 

Q:[I ~e 2 cos(<J ) 2 — 4 >j)] = Q 3 [l -h e 3 cos 4 > 2 ] 

(31) 

Q 2 e 2 sin(<J > 2 — <j>, ) = Q 3 e 3 sin<J > 2 

(32) 

— Q 2 e 2 sin(<J > 2 — <t>i ) = Q 3 e.i sin ef > 2 

(32) 

Dividing equation (32) by equation (31) yields 


Dividing equation (32) by equation (31) yields 


e 2 sin(4»’ — 4>! ) e, sin <J>, 

(33) 

-e 2 sin(c() 2 -((>,) e 3 sin 4> 2 

(33) 

1 + e 2 eos(cf ) 2 — <J>| ) | + e 3 cos <i> 2 

1 - e 2 cos( <t > 2 - <|> , ) 1 + e , cos <|> 2 

Solving for e 2 yields 


Solving for e 2 yields 


e^ sin 

^ — I 

(34) 

— e 3 sin cb 2 

(34) 

s i n ( ct> ^ — <4> i ) — e 3 sin , 

sin(<t > 2 - 4>i ) - e 3 sin cf) i 

Then 


Then 


s i n ( cf> 2 — <(>, ) 4- e d sin ()>-» — sin <()|) 

I +e 2 — — 1 

(35) 

( sin(<J > 2 — <J) i ) + e 3 (sin d) 2 -sin c)> | ) 

(35) 

sin(cj > 2 — cJ>j ) - e 3 sin 4> 1 

sin(cj > 2 — <!>,) - e 3 sin cj). 

From equation (30) 


From equation (30) 


0 _ T sin(c (> 2 - 4 >, ) - e 3 sin cj), 1 


^ ^ T sin (<|> 2 — <*»,) -e 3 sin <j>, ] 

(36) 

[^sin (()> 2 — <J), ) + e 3 (sin d> 2 — sin 4>,)J 

t JDJ 

V2 “ s2l| — -1 

|^sin (<}) 2 — cj), ) + e 3 (sin 4> 2 — sin <t>i)J 

Inserting e 2 and Q 2 into (32) and clearing gives 


Inserting e 2 and Q 2 into (32) and clearing gives 


q sin( 4>2 — <J> i ) — sin <J) 2 = -sin cf> , 

(37) 

q sin((f ) 2 — (f), )-sin $ 2 = -sin cf>, 

(37) 
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where 


Qi -Q3 

Q3 e 3 


( 38 ) 


From (either) equation (37) 

( I — qcos 4>, ) sin + (q sin 4>i)cos 4> 2 = sin cj), (39) 

Notice that equation (39) is satisfied by (Jq = <i> 2 , a solution specifically excluded by equation (36). 
Setting 


I — q cos cf> t 

tan 10 = : — ; 

q sin cpi 


(40) 


so that 


q sin 4>i 

cos uj = — — 

Vq 2 - 2q cos cf> ] + I 


1 — q cos (|>i 

sin to = — — 

Vq 2 — 2q cos cj>| + 1 


allows equation (39) to be written as 


(41) 


(42) 


, • • , sin 4>, 

cos co cos q >2 + sin to sin (p 2 = — : 

Vq 2 — 2q cos 4>| + 1 


or 


cos(cj) 2 - w) 


sin (J)| 

Vq 2 — 2q cos 4>i + 1 


(43) 
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Since, lor any angle 0, 


cos 0 = cos(2tt — 0) 

equation (43) has the two solutions 

, sin 4> | 

<J >2 = co + cos — 

Vq 2 — 2q cos + 1 


( 44 ) 


or 


< 1>2 = 2tt + CO -cos' 


sin 4>| 

Vq 2 -2q cos cj) t 4- 1 


From equation (44), using equations (41) and (42) comes 


sin <}>2 = 


(1 - q cos 4>i) sin <J>| + q sin 4>i V(q — cos 4>1 ) 2 
q 2 — 2q cos 4>i + 1 


(45) 


(46) 


and 

q sin 2 4>, - (1 — q cos cj),) V(q — cos <b|) 2 

COS CD') — ^ ; 

q~ — 2q cos 4*1 4- 1 

Suppose q > cos 4>i, then equations (46) and (47) give 


(47) 


sin <f >2 = sin 

(48) 

cos <j >2 = cos 4>j 

(49) 

But <j> 2 = 4>i is specifically excluded by equation (36). 


Suppose that q < cos <J>|, equations (46) and (47) give 


( 1 — q 2 ) sin 4>, 

sin 4>2 = , 

q~- 2q cos cp j + 1 

(50) 

, 2q - (1 +q 2 ) cos (J>, 

COS 4)2 - , - , , . 

q* - 2q cos 4>i + ‘ 

(51) 
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The conclusion is, then, that if 


q < cos c(>i 


then 


(J>2 


to + cos 


sin 4>i 

Vq 2 — 2q cos 4»i + 1 


From equation (45) 


sin 4>2 


(1 — q c os 4>|) sin 4> i — q sin 4>| V(q — cos 4>i)~ 
q 2 — 2q cos i + 1 


(52) 


COS <t>2 


q sin 2 <J>, + (1 — q cos 4>|) V (q cos cj>i) 2 
q 2 — 2q cos <j)| + 1 


(53) 


Suppose q > cos <t>i, then 

(1 - q 2 ) sin <J>i 

sin 4> 2 = 2 — t I 7 i 

q — 2q cos cp| + 1 


and 


2q - (1 +q 2 ) cos <J>, 

cos 4> 2 = — " — ^ 7 ~~T 

q~ — 2q cos cpi + 1 

Now, suppose cos 4>i > q, then 

sin <f >2 = sin <j)| 


(54) 


(55) 


(56) 
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and 


COS <J>2 = COS <t> , 


(57) 


And the conclusion is that if 


cos <f>| > q 


then 


<|> 2 — 2 "IT + 00 — COS' 


sin 4>| 


Vq- — 2q cos cj), + 1 


since <j>| — cj) 2 is not allowable. 

mm o " r “ alIowab ! e ” is und erstood in the context that <J>, = <J> 2 introduces a singularity in equation 

36) so that Q 2 diverges to infinity. This is, physically, allowed because it corresponds to a parabolic trans- 
ter. I he case ot cj) 2 — 4>j is rejected, however. 

Equations (50) and (5 1 ) [or (54) and (55)] are always valid. Using these, the term sin(4^-4>, ) which 
appears in equations (34) and (36) can be written as 


sin(cj). -<{>,) 2 sin 4>, (cos<|>,-q) 

q--2q cos cf> , + I 


(58) 


so that either ot equations (34) can be combined into 




± e 3 (l -q 2 ) 

2(cos 4» [ — q) — e 3 (q 2 — 2q cos <J>, + 1) 


(59) 


Using equation (37) to define sin <J> 2 and equation (58) to define 
becomes 


sin (<j> 2 -c]>i) equation (36), for Q 2 , 


Q: = 


( Q ' .) 1 

, e (~< 4 2 — 2q cos <J>, + ll j 

\ 1 + e 3 q / l 

3 L 2(cos 4>i — q) J ( 


( 60 ) 
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However, using the definition of q from equation (38) reduces this to 


Q 2 — Q 3 


e 3 (q 2 ~ 2q cos 4>i + 1 ) 
2(cos <f)| — q) 


( 61 ) 


The sign ambiguity present in equation (59) can now be resolved as follows. If the initial and final 
orbits intersect then the perigee of the final orbit must be less than the radius of the first or the apogee of the 
final orbit must be greater than the radius of the first. In the first case 


a, > a 3 (l-e 3 ) 


or 


1 1 
> — 

a 3 (l-e 3 ) a. 


or 


1 +e 3 > __1 

a 3 (l-e 3 2 ) a, 


or 


Q3(l+ e 3) > Ql 


so that 


Q3 e 3 > Qi - Q 3 


or 


1 > q 

In the second case, 

a 3 (l +e 3 ) > a, 
which leads to 
q > -1 
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Thus, the criterion for intersecting orbits is that 


-1 < q < 1 (62) 

The numerator of e 2 is, then, certainly positive. Consider the denominator of e 2 in equation (59). 
Temporarily defining 

f(4>,) = 2(cos 4), -q)-e 3 (q 2 -2qcos 4> , + 1) (63) 


the extrema of f occur at 0 or 7T so that 

f(0) = ( 1 q) [2 - e 3 ( 1 - q)] > 0 
f(Ti) = -(1 + q) L2 + e 3 (l+q)] < 0 


af 

a<})| 


—2 sin <±» | ( 1 + e 3 q) 


For |q| < 1 , e 3 < I , 


af 

a^i 

af 


o 


o 


(0 =£ 4> t «£ TT) 


(7T ^ 4>| s; 27T) 


(64) 


(65) 


Also, f will vanish (i.e., e 2 will experience a singularity) whenever 


cos 4>| 


2q + e 3 (l + q~) 
2(1 + qe 3 ) 


i.e.. 


. -i f 2q e 3 ( 1 +q 2 )l 

= C ° S l 2( 1 + qe V) J 


( 66 ) 


(67) 
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or 


<J>i = 2tt - cos 


_i | 2q + e3(l + q") 
L 2(1 +qe 3 ) _ 


( 68 ) 


The singularities of Q 2 and e 2 are not coincident. They may be ordered as follows: 


q 2 < 1 =► q 2 e 3 < e 3 =► 2q 2 e 3 < e 3 + q 2 e 3 — e 3 ( 1 + q 2 ) 


Thus 


2q + 2q 2 e 3 < 2q + e 3 ( 1 + q 2 ) -*=► 2q( 1 + qe 3 ) < 2q + e 3 ( 1 + q 2 ) 


so that 


2q + e 3 (1 +q 2 ) (69) 

2(1 +q e 3 ) 


Thus, if 4>! begins at zero and progresses through increasing values toward -it, and if q > 0 the first 
singularity will be encountered at a value given by equation (65) and, subsequently, the singularity cf>j = 
cos' 1 q will be encountered. (In physical terms, e 2 will diverge before Q 2 diverges. ) If we proceed from 2 tt 
through negative values toward t r, the same sequence occurs. 

Let the upper half of the polar orbital plane (0 cf> | =3 it) be divided into three regions, labeled I, II, 
and III, defined as follows: 


I. 


II. 


III. 


4) I < cos 


.1 ~ 2q + e 3 (1 + q 2 ) 
. 2(1 +e 3 q) 


4) | > cos 


1 2q + e 3 (1 +q z ) 


2( 1 + e 3 q) 


4>| > cos 


_] r 2q + e 3 (1 + q 2 ) 
2(1 +e 3 q) 


and 4 > i cos 1 


and 4>, < cos' 1 q 


and 4>, > cos' 1 q 
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The goal is to define physically realizable orbits over as many of these regions as the constraints allow. 
From f(0) > 0 and equations (64) and (67), e 2 in region I must be positive, so choose the sign in equation 
(59) such that 


Ct — 


+ e 3 ( 1 - q 2 ) 

2(cos cj), - q) - e 3 (q" - 2q cos <j>i + 1) 


and p = — <j>|. Also, for |q| < i,Q 2 (0) > 0. Inserting (66) into (61) gives 


q 2 


| cos" 1 j” 


2q + e 3 ( 1 + q 2 )' 


2(l+qe 3 ) 


= 0 


(70) 


(71) 


(i.e., parabolic transfer). Thus, region 1 is an area of valid orbital transfers, and similarly for the mirror 
image in the lower half of the plane. 

In region II, Q 2 decreases from 0 to negative infinity, always staying negative. This region, then, 
does not correspond to physically realizable transfers since Q 2 must be positive. 

Region III again finds Q 2 to be positive since 


Q 2 (tt) > 0 


and 


dQ 2 

dQ<j> 


< 0 


(0 <})| s£ tt ; <|>, cos" 1 q) 


Since Q 2 undergoes an infinite discontinuity at (J>] — cos" 1 q, these results are consistent. 

The only remaining question is the behavior of e 2 in region III. From (59), certainly 

q 2 - 2q cos <j>| + 1 > 0 

and since <j>, > cos" 1 q. 


cos <f>| - q < 0 
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so that 


2(cos 4>| - q) - e 3 (q 2 - 2q cos cf> j + 1) < 0 


Since I — q 2 > 0, the conclusion is that e 2 will be positive if 


- e 3 ( 1 - q 2 ) 

2(cos <J)| - q) - e 3 (q 2 - 2q cos cf> j -(- 1 ) 


( 72 ) 


from which the second of equations (34) shows that p = tt - <f>. 

The information gathered on the regions can now be summarized pictorial ly as follows: 


e 2 from eq. (70) 
Q 2 from eq. (61 ) 
P = - <t>i 

Forbidden zone 
No physically 
realizable 

tangential transfers 

e 2 from eq. (72) 
Q 2 from eq (61) 
p = TT - 4»i 

♦ ♦ 4 


<J>, = 0 


e 3 (1 +q) 

2 - e 3 ( 1 - q) 


Qi = Q> [i - C3(l 2 — ] 


c|) | = cos' 1 

f 2q + e 3 ( 1 4- q 2 )"l 

L 2( 1 + q e,) J 

e2 = x 


Q 2 = 0 



<j)i = cos' q 


e-, = I 


0 - = 


<1>I == IT 


e-> = 


e 3 (l -q) 


2 + e 3 (1 +q) 

Q: * Q,[« + ^]| 


These are valid for any intersecting orbits with ei =0 and e 3 ^0, but no discussion has yet been 
offered for non-intersecting orbits (|q| > 1). That case is considerably easier to deal with than the intersect- 
ing case. The equation for e 2 , equation (59), is valid but now, since |q| > 1 , there can be no singularity in 
e 2 . This is because equation (66) cannot be fulfilled (|cos <{>i| *£ 1). 
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From the bounds on f, following equation (63), 


(1 - q) (2 - e 3 ( I - q)] 3* 2(cos 4>, - q) - e 3 (q 2 - 2q cos <J>| + 1 ) > - (1 + q)[2 + e 3 (l +q)J 

If q > 1, then f < 0. If q < -I then 1 > 0. This yields a simple rule since 1 - q 2 < 0. 

If q > 1 , choose 

e 3 ( 1 - q 2 ) 

e-> = , 

2(cos 4>i — q) — e 3 (q- - 2q cos (J>, + 1) 


if q < — 1 , choose 


- e 3 (l - q 2 ) 

2(cos | - q) - e 3 (q 2 - 2q cos 4> t + 1) 


The calculations on Q 2 are more delicate. From equation (61), if q > 1 it is obvious that Q 2 > 0 for 
all <J) 2 and e 3 . It q < - 1 then additional considerations are needed. From equations (38) and (II) comes the 
relationship 


Qi _ Q3 
q Q3 e 3 


1 a 3 (l-e 3 2 ) _ a 3 (l - e 3 2 ) - ai 

£3 e 3 

a 3 (l - e 3 2 ) 


(73) 


so that 


3l 


a 3 ( 1 - e 3 2 ) 
1 +e 3 q 


(74) 


or 


a 3 


a,(l +e 3 q) 
1 - e 3 2 


(75) 
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In each case both a | and a 3 must be positive along with 0 =£ e 3 < I . But combinations can be chosen 
which violate this physical constraint. (There is no treatment here of the initial or final orbits being 
parabolas or hyperbolae.) The requirement, then, is that 


1 + e 3 q > 0 


or 


— 1 >q> — l/e 3 , 0 < e 3 < I 

Let e be a real positive number defined according to 

q = e 2 - l/e 3 . 

From the extremal values on Q 2 comes the inequality 

Q , s Qj « o, [, + Mi±S>] 


or 


(Q 3 /2) 11 -e 3 (l -e 2 )] ^ Q 2 < (Q 3 /2) [1 + e 3 (l +e 2 )] 


so that Q 2 will be positive. 

The conclusion that is reached is that, for the present case, there are no restrictions on tangential 
orbit transfers if the orbits do not intersect. 
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Third Singular Case: Initial Eccentricity # 0, Final Eccentricity = 0 

This case is not identical to the previously treated case. The analogues of equations (24) through 
(27) are 

Qi( I + e, cos <t>,) = Q : [l + e 2 cos(4>, + p)] (77) 

Q, e, sin <J), = Q 2 e 2 sin(<J>i +p) (78) 

Q 2 | 1 + e 2 cos(4> 2 + p)] = Q3 (79) 

Q 2 e 2 sin(4> 2 + p) = 0 
from which comes, immediately. 


p = -ch (81) 

p = 7T - <J) 2 . (8 2) 

In a manner exactly parallel to prior developmerits the following equations are obtained. 

If P = — (}> 2 


ei sin <J>, 

sin(4> | - <|> 2 ) - ei sin ct> 2 


(83) 


If P — TT - 4> 2 


-ej sin (j>| 

sin( 4> 1 - (}) 2 ) - ei sin c() 2 


(84) 
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In either case 


f sin(4>i - 4> 2 ) ~ 
^ 2 sin(4> j - 4> 2 ) + £i ( 


e t sin 4); J 

sin 4>i — sin 4 * 2 ) J 


The definition of q is changed, tor this case, to read 


q 


Qj - Qi 
Qi e i 


(85) 


( 86 ) 


then 

(1 — q cos 4*i) 4*2 t (q sin 4 m) c °s <J>2 = s * n 4*1 

The value of to from equation (40) has the same functional form so that equations (44) and (45) are still valid 
for <J> 2 ; similarly, for the values of sin <J> 2 and cos «J> 2 as given by, say, equations (50) and (51). 

Since the functional form of e 2 differs, the former equation cannot be expected to hold true. Now 


2 sin 4>i (q - cos <l>i ) 
sin(<J,, - <M = q- _ 2q cos + ] 


(87) 


so that from equations (83) and (84) 


±e, (q 2 - 2q cos ! + P (88) 

2(q - cos 4>,) — ei (1 - q 2 ) 


The equation for Q 2 , (85), becomes 

T e,(l - q 2 ) 1 

Q 2 - Qi 2 ( q_ cos <*>,) J 


(89) 


Next, the sign of e 2 must be determined. This will be done in the same manner as before, first for the 
case of intersecting orbits. 
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Assume that the perigee of the initial orbit is less than the radius of the final orbit 


ai ( 1 — e | ) < a 3 


_ 1 _ I 

a 3 a , ( 1 — ei) 

_ L< 1 + e i 

a 3 a,( 1 — e, 2 ) 

Q> < Q, (1 +e,) 


Therefore 

q < 1 

Next, assume that the apogee of the initial orbit is greater than the radius of the final orbit so that 


ai( 1 H- e | ) > a 3 

J_ 1 

a 3 > a, ( 1 +ei) 


1 

> 

a 3 


1-e, 

aid ~e, 2 ) 


Q 3 >Q, (1 - e,) 


Qa-Qi 
Qi e. 


> -1 


so that the familiar bound is still 


-I <q< 1 


( 90 ) 


22 



From equation (88), e 2 will experience local minima at 4>i = OandcJ>i — tt. Furthermore, e 2 has a 
singularity at 


4> ! — COS 


-1 


r e , ( 1 - q 2 ) 1 

L q — ~J 


At <J>, = 0, 


(91) 


= ±e,(q - 1) 

2 + C|( 1 + q) 

So to ensure a positive e 2 , the negative sign is chosen. The same sign choice continues until e 2 diverges. At 

4>i = 77 


__ ±ei(q+ 1) 
62 2 + e,(q - 1) 


(92) 


so that in the locality of tt, the positive sign is the proper choice. 
Q 2 [eq. (89)] will experience a singularity at 


4> | = cos' 1 q 


The singularities may be ordered as follows: 


1 - q 2 > 0 


ei 

2( I - q 2 ) 


0 > 


-ei 

2(1 -q 2 ) 


q > q 


£i 

2(1 -q 2 ) 
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Thus, if <J>, begins at 0 and increases toward tt, the first singularity that will be encountered occurs 
in Q 2 . At that point Q 2 diverges and becomes negative. It continues as negative until the singularity in e 2 is 
encountered. Thus, the “forbidden” region is determined by 


cos' 1 q cj), 


cos 


q 


q- 


£i 

2(1 -q 2 ) 


(93) 



As before, there are no restrictions on the value of <t> , for transfer if the orbits do not intersect. 

The General Case: Initial Eccentricity £ 0, Final Eccentricity ^ 0 

The prior three sections assumed that either or both of the eccentricities of the initial and final orbits 
vanished and this yielded much simplification. The general case will now be treated. Equations (12) 
through (15) constitute the natural starting point for the development. 

Dividing equation (13) by equation (15) yields 


sin(<j> | + p) _ Q, ei sin 4>| 

sin(<t> 2 + p) Q, e 3 sin(<J> 2 + 0) (94 ) 
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Defining 


Qi e, 

7 = q7^7 

allows equation (94) to be written as 

sin(4>, + p)sin(<j> 2 + ft) = 7 sin 4>i sin(4> 2 + P) 


Expanding the functions which involve p. gathering, and solving for tan p gives 


tan p 


y sin <\>\ sin d>: - sin (j>i sin(d >2 + ti) __ 
cos ((), sin(<h + fl) - 7 sin cos ^ 


sin d> 




<7 


$) lan d>: - sir> ^ 


c() t cos 


ft tan <J) 2 + cos <f>, sin ft - y sin <t>. 


From equations (12) and (13) comes 


e 2 sin( 4>| + p) _ e i s ‘ n ^i 
1 +e 2 cos)^! + p) l+ei cos cj>i 


Solving for e 2 


ei sin 4q 

e2 sin(4>i + p) + C| sin p 


Multiplying equation (12) by sin(<|>,+p). equation (13) by cos(«t>, + p), and adding 


q 2 = 


Q | [ sin(cf) t + p) + e, sin p] 
sin(4> i + p) 


Equations (98) and (99) can now be used to eliminate e 2 and Q 2 from equation (14). 

Q ^e! sin p + sin(4)| + p)l l~. e i sin 4>i cos(4> 2 + P]_~l _ n,[ I + p-,cos((b-. + ft)l 
sin(4> ! + p) L sin(<J>, + p) + e, sin pj 


(95) 

(96) 

(97) 


(98) 

gives 

(99) 
This yields 

( 100 ) 
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Clearing shows that either 


e, sinp + sin(4>, +p) = o 


( 101 ) 


or 


Qi{sin(d>, +p) + e (l sin p + sin 4>, cos(<(m + p)|} = Q,[| +e,cos«|> ! + a)lsin«|>, +p) ( | 02 ) 

Comparison berween equations ( 100) and (98) leaves only equation (102) as a valid alternative. 

Now using the identity 

sin(<f>, -t- p) = cos pfsin <J>, + cos , tan pj (|03) 

and substituting for tan p from equation (97) into equation (103) gives 

sin(t|) l + p) =. 7 sin 4), sin((fj 2 - 4>,) cos p 

cos <J) 2 | cos 4>| cos O tan <J> 2 + cos 4> ( sin -0 — -y sin c(>|] (104) 

Similarly, since 

cos(4> 2 + p) = COS p [cos c|) 2 - sin (f> 2 tan p] 

= t~ 7 s * n foi + cos(([)| ~ <{>2) sin(<j) 2 + ft)] cos p 

cos <|) 2 (cos (f> 1 cos p tan (J> 2 + cos <p, sin - y sin <j> r ) (l05) 

A portion of the second member of equation (102) can now be constructed as 

sin p + sin <J>, cos(<t> 2 + p) = cos p {tan p + Z 1 Sin ^ + cos (4>i - <t>2> sin(cj) 2 + ft) ^ 

1 c °s 4*2 (cos 4>, cos tan <|) 2 + cos <f>, sin d -7 sin 4>,)/ 

(106) 

Again substituting for tan d from equation (97) and combining fractions gives 
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( 107 ) 


sin p + sin 4>i cos(4>2 + p) 


cos p sin 4*i ( 7 (sin <j> 2 — sin <|>i) — sin(4> 2 + ft)[ 1 - cos(4>i - 4) 2 )]} 
cos c{) 2 (cos 4>, cos ft tan 4>2 + cos <f> j sin ft - 7 sin 4>i) 


Next, substituting from equations (103) and (107) into equation (102) and cross-multiplying by 
cosc|> 2 (cos | cos 0 tan 4>2 4- cos 4>i sin ft - 7 sin cj>i) gives 


Qi ei cos p sin {7(sin 4>2 - sin 4>i) - sin(4>2 + ft)[ 1 - cos(<J>, - 4> 3 )] + 

7 sin (j), sin(4> 2 — 4>i) cos p 

cos 4 > 2 [cos (J)| cos ft tan 4>2 + cos 4>i sin 0 - 7 sin 4> 1 ) 

(108) 


= Q 3 1 1 + e 3 008(4)2 + 13)] 7 sin 4h sin(4> 2 - 4>i) cos p 


Require that none of cos p, sin 4>i , or cos 4> 2 vanish. Then substitute for 7 from equation (95) on the right 
side of equation (108), then expand equation (108) to obtain 


e, e 3 {7(sin 4> 2 - sin 4>,) - sin(4> 2 + ft) [ 1 - cos(4>, - 4> 2 )]} + e 3 7 sin(4> 2 - 4+) 

(109) 

= e! sin( 4>2 — 4> 1 ) + £i e 3 cos(4>2 + ft) sin(4> 2 - 4>i) 


Eliminate the remaining 7 terms from equation (109), expand the multi-angle expressions and gather on 
sin 4>2> cos 4*2 to obtain 


[Qi e i - Qi e 3 cos ft + (Q| - Q 3 ) cos 4+] sin 4>2 - [Q3 e 3 sin ft + (Q, - Q 3 ) sin 4>i ] cos 4*2 
= Qi e| sin 4> 1 - Qi e 3 sin(4)| + ft) 


( 110 ) 


For convenience, define 


Pi = Qi C|-Q 3 e 3 cosft + (Qi-Q 3 )cos 4>i = Q 3 e 3 (7-cosft + qcos 4>i) (111) 

p 2 = Q 3 e 3 sin ft + (Q|-Q 3 )sin 4>i = Q 3 e 3 [sin ft + qsin 4+] (1 12) 


* Letting sin 4>i = 0 can provide a valid transfer. 
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p 3 = Qi e i sin <t>|-Q 3 e., sin(cj>, + f>) 


(1 13) 


where q is given by equation (38), not by equation (86). 

Equations (111) through (113) are not independent as can be seen from 

p, = Pi sin <t>i -P 2 COS 4>i 

so that equation (110) can be written as 

P! sin 4*2 — P 2 cos 4>2 = pi sin <J>i -p 2 cos 4>| 

Dividing equation (1 15) by Vp, 2 + p 2 2 and defining a new to by 
to = tan" l (pi/p 2 ) 

allows equation (115) to be written as 

cos t |>2 cos to - sin 4>2 sin co = cos 4>| cos 00 — sin <|>i sin to 
or 

COS(tJ >2 + to) = COS(4)| + co) 

Equation (117) has two solutions. The trivial solution is 
<j) 2 = 4>i 

but the desired solution comes via the identity 
COS(c|) 2 + to) — COS(2tT - t )>2 - to) 


(114) 


(115) 


( 116 ) 


(117) 
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Then 


2tt — 4>i — <*J — 4>i + w 


or 

cj) 2 = 2 (tt - a)) — 4>i = 2 [tt — tan ‘(p|/p2)J — 4*i 


(118) 


Equation (118) is the generalized form of either equations (44) and (45). 

The next step is to utilize equation (118) to eliminate <J> 2 from equation (97) in order to obtain an 
expression for p which depends only upon <t>i * 


Defining the variables 


a = tan 


sin -9 


cost) 


(119) 


r = Vy 2 — 2y cos f) + 1 


( 120 ) 


naturally leads to 


sinf) = r sin a 

(121) 

y — cos 0 = r cos a 

(122) 

Then equation (97) may be expanded as 



ry^hor wishes to thank J . R Redus of Marshall Space Flight Center, Huntsville, Alabama, for point- 
ing out the following transformation. 
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tan p = sin <f>i 


( 123 ) 


7 sin (j) 2 - sin 4> 2 cos ft - cos <j> 2 sin {) 

cos 4>| sin 4>2 cos O + cos 4>i cos <J> 2 sin 0 - 7 sin <|>i cos <b 2 

r sin 4> : sin(4> 2 - a) 

7 sin(<() 2 - <j>|) - r cos(cJ) 2 - a) 


Now, substituting for cj> 2 from equation ( 118 ) yields 


tan p 


rsinfoi sin [2-rr — 2fa>- <j), - a] 

7 sin[ 27 T - 2a> — 24>, ] - r cos <f> , sin[ 27 T- 2 io-<t>i ~aj 


_ -r si n (j), sin(2a) + cf), +a) 

-7 sin[2(co + c}>| )] + r cos <f>, sin(2oo + <J>| +a) 

For convenience, define 


R = Vr + 2 rqcos( 4 >| - a) + q 2 


Then 


r sin a + q sin 4), 

cos a) = - 

R 


r cos a + q cos cb. 

sin u> = 

R 


so that 


COS(tO + <j) | ) 


- r sin(ct>] - a) 

R 


( 125 ) 
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r cos(<}>i - a) + q 
sin(co + 4»i) = d 


cos(co + a) 


q sin(4> i - a) 

R 


sin(a) + a) 


r + q cos(c{)| - a) 

R 


Thus, 


(q 2 - r 2 ) sin(4>i - a) (126) 

sin(2w + <|> , + a) = sin(co + 4> i ) cos(w + a) + cos(w + 4> , ) s 111 ^ + a) R 2 


_ -2r sin(4>i - «) [r cos (4>) - a) + q] 
sin[2(<i) + 4>i)] = 2 sin(u) + <J>i)cos(a) + 4>i) - R 2 


(127) 


.. c m? 61 and (127) into equation (124) gives 
Inserting equations (126) ana t 

. ,, v /p 2 -i r - r(q 2 - r 2 ) sin 4> t 

lan P = [zn(r cos«h - “> + c l' + r COS (q2 " r ^ 


(128) 


S arctangent function is ntu.ti-valued i, wili no, yield the sante value for 


tan' 1 


f(x) 

g(x) 


as it does for 


tan 


- f(x) 

- g(x) 
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Equation (128) can safely be written as 


tan p 


— r (q 2 — r 2 ) sin 4q sign( sin( 4? | - a)] 

{2ry[r cos(<J), - a) + q) + r cos 4>i(q' - r)} sign[sin(4», - a)] 


( 129 ) 


Eliminating a and r from equation (129) gives the final form of equation (129) as 


tan p 


(7 2 - 2~y cos fl + 1 - q 2 ) sin <j>, signfsin(4>, - a)l 

[( 7 ~ + q- - 1) cos <J>| + 2y (q + sin fl sin 4>,)] sign[sin(<|>, - a)] 


(130) 


= iX_ 2y cos 0 + | - q 2 ) sin 4>, sign[y sin 4>, - cos(<j), - ft)] 

[(T + q " ~ !) cos + 2 ^(q + sin ft sin 4>,)J signfy sin 4>, - cos(4>, - ft)] 


The machinery is at hand to obtain an expression for Q-, in 
equation (99) comes 


terms of fundamental variables. 


From 


Q : = q, i + 

t Sln *i COS p + cos <J>, sin p 

= Q, It + £i lan £ * 

( sin 4>i + cos <f>| tan p j 


( 131 ) 


Using equation ( 1 30) (and 
expression 


ignoring the sign terms which make no difference here) yields the 


sin 4>, + cos 4), tan p = 2-ysin <J>, L 7 cos 4>i + q - cos(<f>, + 0) 1 

L(y 2 + q 2 - I) cos cf), + 2y(q + sin 0 sin 4),) J (132i 

Then 

Q, = Q,jl + e|( T ~ 2y cos 0 + I - q 2 ) ) 

I 2y [y cos 4>i + q - cos(4>, + 0)] j (133) 

is a fundamental expression for Q,. 
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The final equation which is needed is the expression for e 2 which involves only fundamental vari- 
ables, i.e., a non-sequential calculation, I have been unable to derive an expression which is aesthetically 
satisfactory, but such an expression is not required in what lollows. 

The “forbidden regions” which constitute the focus of this paper exist in the general case and can be 
analytically defined by enforcing the following conic section restrictions: 


|(0 < e 2 < 1 ) and (a 2 > 0)] 
or [(e 2 > 1) and (a 2 < 0)1 
or [(e 2 = 0) and (a 2 > 0)1 

or |(e 2 = I) and (a 2 = oc )] 


( 134 ) 


It will be shown that e 2 and a 2 can take on any value on the extended reals and the regions which violate 
equation ( 134 ) give rise to forbidden regions. 

As in the special cases, the apparent dependency of sign(e 2 ) on sign[sin(<J>i )] is illusory [equation 
( 98 )]. Consider these points at which e 2 will experience a singularity, i.e., the solutions of 


sin(4>| + p) + ei sin p = 0 


035 ) 


which expands to give 

sin(<t>i + p) + C| sin p = [(7" + q~ — 1 ) cos 4>i + 2y(q + sin 0 sin 4>i)] sin 4 m 

+ [(y 2 -2ycosf) + I -q 2 )]cos cj), sin 4> 1 + e i (y 2 -2ycosf) + 1 -q - )sin 4 >i = 0 

Since sin <J>, is assumed not to vanish, gather on sin <J>|, cos <J>,,to yield 

„ . , i ei[2y cos ft + q 2 - y 2 - 1] -2yq ( 

(y-cosf))cos c()| + sin 0 sin cpi = ) l 


( 136 ) 


( 137 ) 


Utilizing the definitions of ol from equation ( 119 ) and defining 


£ = cos' 


e 1 [2-7 cos ■d + q 2 - 7 2 - 1] - 2yq 
2y V 7 2 - 2y cos 0 + I 


( 138 ) 
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then equation (137) has two solutions 


<J>! = OL + £ 

4>i = 2tt + a - £ 


(139) 


as singular points of e 2 . 

Another complicating factor to the orbit transfer arises when the numerator and denominator of 
equation (130) vanish simultaneously. These points can be isolated as follows. 

Assume that the numerator of equation (97) vanishes [this is more convenient than working with 
equation (130)]. Then 


(y-cosO)tan <J) 2 -sin 0 = 0 


(140) 


so that 


tan c]> 2 = tan a ( 141 ) 

If the denominator of equation (97) also vanishes 

cos 4> | cos 0 tan cj) 2 + cos 4> | sin 0 - y sin <j> ] = 0 (142) 

Eliminating <f> 2 between equations (141) and (142) gives 
. sin 0 

tan <j>, = -= tana (143) 

y - cos O 

Since 


tan 4>| = tan(7T + <t>|) 
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the two solutions of equation (143) are 


4>, = 4>2 = it + a 

which are also the points at which sin(4>i - a) = 0. Thus, when the numerator and denominator ot equation 
(97) [or (130)] vanish, sign[sin(4>, - a)] loses meaning since zero is unsigned. 

Both the numerator and denominator are continuous through the points given by equation ( 1 44) and 
the angle p will experience a jump discontinuity of it. The behavior of e 2 at this point is demonstrated as 
follows. Let <JT and <J> + represent 4> to the left and right of the solution given by equation ( 144). Then 


e 2 (<|> ) 


ei sin cj> 

sin(4r + p) + e, sin p 


e 2 (<}> f ) 


e| sin c}> 4 

sin(4> + + p ± tt) + e, sin(p ± tt) 


sin(4»+ p ± tt) + e, sin(p±7r) = cos(±tt) [sin(cj) + p) + e, sin p) 
= - [sin(4> + p) + e, sin p] 


|e 2 (4> )| = - |e 2 (4> + )| 


This far, e 2 has exhibited two unbounded discontinuities [see equations (139)] and two bounded discon- 
tinuities [see equations (144)]. Additional information on e 2 comes from the observation that it 


e 


(146) 


(i.e., parabolic transfer) and it 


4>, + p - 0 


( 147 ) 
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then 


e-. = — I 


( 148 ) 


To isolate the value of <t> | which corresponds to equation (145). observe that, necessarily 


tan p = - tan 


Inserting this result into equation (130) gives 


(y 2 + q 2 - 1 ) cos <f)| + 2y(q + sin Osin <j>|) = -(y 2 -2ycosO + 1 -q 2 )cos <j). 


( 149 ) 


Gathering, utilizing a once more and defining 


Vy 2 - 2y cos 0 + I 
then equation ( 149) has the two solutions 
<f>i = a + 4 

<j>i = 2-ir + a - £ 


(150) 


(151) 


The physical identification of the angles given by equation (151 ) is worthwhile. Suppose that equa- 
tion (118) is valid and. simultaneously, the trivial solution is valid so that 

<l>i = <t>2 


This corresponds, of course,* to orbital intersection. Then equation (1 18) requires that 
tan <t>, = - pi/p 2 


* The same results can be obtained by equating the left member of equation (12) to the right member of 
equation (14). 
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Inserting p, and p 2 from equations (111) and (112) yields 


. , 053) 

(-y- cos 0) cos (J), + sin 0 sin <pi — q 


which is the reduced form of equation (149). Thus, at an orbital intersection, e : - +1- 

The orbital transfer plane is cleaved into two ^ ”„ S “uch points have 

it is reasonable to suspect that there may be four P aces w e I 2 ' , ions ( 1 2) and ( 1 3) yields the 

been isolated. Two others can be obtained as follows. Setting M - I m equations ( 12) ana y 

pair 


Q, (l+e, cos 4>i) = Q 2 M ± cos(<j>, + p)l 
Qi ei sin <J> i = — Q: sin(cj>| -H p) 


Isolating sin(4>i + p) 


and cos(<J>, +p). then squaring and adding the resulting equations gives 


Q , ( 1 + 2e i cos 4), + ef) = 2Q 2 (1 +ei cos 4>i) 


(154) 


Inserting Q 2 from equation (133) and cancelling like terms gives 


fe, 2 (q 2 + y cos # - 1) - y(y ~ cos ft)] cos 4>, + y(e, 2 
= Ci('y 2 -2'y cost) + 1 -q") _ q7( e r~ H] 


1) sin f) sin 4>i 


(155) 


Defining 


a 


, y(e, 2 — 1 ) sin -ft 

tan e, 2 (q 2 + y cos d - 1 ) - y (y - cos d) 


(156) 


and 


t = cos' 


Cl ( y 2 - 2~y cos ft + 1 — q ) ^ — 

V(e, 2 (q 2 + "Y cos 0 - 1 ) - - cos ti)| 2 + [y(e, 2 - D sin ()F 


(157) 
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then the other two points of unity for e 2 are given by the solution of equation (155) as 
4> i = cr + t 

4>i = 2tt 4 - a - t * 


•men.' W Th' n th,S SeCt,0n the behavior of e 2^i> has been a focal point while Q,(<j>,) has received less 
attention. The next step is to examine Q 2 (cf>,) at each point which has proven to be of interest in e 2 «K). 


The first such point pair is given by equations (139). The eas 
writing equation (137) as 


easiest way to make the substitution is by 


jMT - 2y cos 0 + | - q 2 ) 

~~ 2y = -q-(7~cos^)cos <J>| -sin 0 sin <f>, 

and substituting the left side of equation (159) into equation (133) which immediately yields 


(159) 


Q 2 = 0 


( 160 ) 


at this point. This value then shows that a hyperbola of infinite eccentricity has resulted, i.e. , a straight line, 
e, From^equah io'™ (T^Tco me s ^ ^ eqUati ° nS ( ,44) * ^ P ° mt ° f fi " ite J Um P discontinuity in 


sin cf> | = 


sin 0 


Vy 2 - 2y cos f) + 1 


cos <j> I = 


± ( y - cos 0) 
Vy 2 - 2y cos -0 + 1 


so that equation (133) becomes 


Q: Qi {1 ± (e ( /2y) (Vy 2 -2ycos O + 1 + q]} 


(161) 


where .he positive sign corresponds to .he 4>, = a and .he negative sign corresponds to = rr + a . 
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The first case where |e 2 | = 1 comes from equations (145) and (147). In either case, equation (99) 
shows a divergence of Q 2 . 

The remaining case occurs when <J), is given by either of equations ( 158). The actual value of Q 2 at 
this point is not difficult to determine, numerically. But a single equation which gives that value is quite 
messy. Since the coefficient of Q 2 in equation (154) cannot vanish, this point is a regular point. 

The following table. Table I, summarizes the information presented here on the behavior ot e 2 (<J>,) 
and QUO- 


TABLE I 


Governing 

Equations 

Value of <{)| 

Value of e 2 (<J>i) 

Value of Q 2 (4>i) 


a + { 

Diverges 

0 

139 

277 4- a - £ 

Diverges 

0 



±ei sin 4 >i 



a 

Qi { 1 4 (e,/ 2 y) [V-y 2 - 2 , y cos i) 41- q|} 


sin(<j)i 4 p) 4 ej sin p 

144 

77 + a 




(Jump Discontinuity) 

Qi { 1 — (e |/2-y) (V 7 2 - 27 cos 0 4 1 4 q]} 


a 4- £ 

± 1 

Diverges 

151 



(orbital intersections) 

277 4 a - ^ 

± 1 

Diverges 


O' 4 T 

± 1 

Regular Point 

158 

i 2 tt 4 a - t 

± 1 

Regular Point 
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It is very convenient to order the angles which are presented in Table I. That is, for a given input of 
the primitive variables, predict the order in which the important points of will occur. It must be realized 
that certain conventions are inherent in the derivation of the prior equations - i.e. , all angles must be chosen 
in the range |0,2tt], the principle value of cos '( ) will lie in the range [0,tt], etc. These conventions follow 
from standard computer implementations of the derived equations and we make use of them in the follow- 
ing discussion. 

Since all angles are calculated in a modulo 2tt system, any angle in column 2 of Table I can exceed 
any other angle. Arbitrarily, but without loss of generality, choose [0 < a < it]. Then 


a + tt > a > 0 


Because of the principle value condition, equation (138) will yield 


tt > £ 2= 0 , 


so that 


2tt > 2-TT — l > IT 


and 


2-TT + a>2Tr + a — ^>TT-l-a , 
since 2tt + a = a and tt + a = cx-tt 

a > 2 tt + a — £>ir + a = ai — tt . 

The next ordering occurs between the angles £; and £. For £, to be real, equation (150) shows that 
— q ^ Vy 2 — 2y cos 0 + 1 . 

For convenience, abbreviate, as before, 
r = Vy 2 — 2y cos 0 + 1 . 
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Then 



so that 


ei(q 2_ r 2 ) 55 0 * 


and 


ei (q 2 -r 2 ) - 2-yq =£ -2^q 


or 


ei(q 2 ~r 2 ) - 2-yq < _ S_ 
2-yr 1 


so that 


cos*'( — q/r) > cos 


e i (q 2 — r 2 ) - 
27r 


From equations (150) and (138) comes 


^ i 


£ + ct =S £ + a 


2'ir + a — £ ^ 2 tt + a — ^ 


The angle ordering, thus far, appears as 


a<a + £<a + £<a + 7T<2-rr + a- £<2TT + a- £ 


By construction, the remaining two important angles c + t and 2 tt + c — t can only occur in one of 
two locations. o- + Tor2-iT + CT-T must occur between 2 + it and 2 tt + a - £ or between a + £ and a + ir. 
Whichever angle is established tor the first slot determines that the second angle will occupy the second 
slot. Thus, we need only establish a single inequality for a chosen configuration to determine a complete 
ordering. Choose a = ir so that ot-TT = 0. From equation (1 19), {> = 0. If 1 16 = 0, then ( 156 ) shows that 
a = 0. Now because of principle value considerations, tt 3= t so that 


2tt 3= 2t 


or 


2 it — T 


T 


and 


2tt + o — t ct + t 


Thus, the cyclic ordering can be written as follows: 


(a + TT < a + T < 2 tt + a — £ < 2 tt + a — £ < a < a + £ < a + £ < 2-rr + a - T) mod 27T 


The computation algorithm is quite simple. Choose the minimum value of any of the important 
points and arrange the rest in accord with the above ordering. 

The above theory has shown that certain regions of intersecting orbits preclude tangential transfers. 
It is intuitively obvious, however, that these forbidden regions cannot occur in the case of non-intersecting 
orbits. It is important to demonstrate that it the orbits do not intersect then there are no discontinuities in 

e 2 (4*i )- 


Consider the following, since Q, and Q 3 are both positive 


Oi > O.-Q3 
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so that 


Qi e, (Qi-Q .0 £i 
Q3 63 Q3 

i.e., 

7 > e, q 

27 > 2e, q = e, q + e, q . 

Now suppose lhat an interseclion does not occur. Then equation ( 150 ) shows that £ will not be real, i.e., 
(-q > r) or (-q < -r) ■ 

Thus, if -q > r, 

27 > e, q - e, r = e, (q-r) 

and 

0 > (r + q) 

so that 

e,(q-r)(q + r) > 27(r + q) 
or 

ei(g 2 -r 2 )- 27 q > { 

2yr 
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Conversely, if — q < -r then q-r > so that 


e i(q~ ~ £~) — 2yq 


2yr 


< - 1 


From (138) it is now apparent that e 2 has no infinite discontinuities in the real plane. It is also easy to 
demonstrate that the t.n.te discontinuities are precluded for non-intersecting orbits. In equation (130) the 
numerator cannot generally vanish via r 2 = q 2 so, if the only way it could vanish is via sin tj>, = 0 If this is 
the case then the denominator will become 


f = ±(T + q 2 - 1) + 2yq 


which can be written as 


f = ± (y ± q) 2 -+- 1 

Suppose that f does vanish. Substituting for y and q in terms of primitive values shows that one of 
the lour conditions 


a i ( i + e i ) = a 3 ( I + e 3 ) 
a i ( I c i ) = a 3 ( 1 + e 3 ) 


a | ( 1 + e i ) = a 3 ( 1 — e 3 ) 
ad 1 e | ) = a 3 ( I -e 3 ) 


which correspond to osculating contacts between the orbits, i.e., single impulse transfers. 

With respect toQ 2 (c|)|), this function will diverge if the denominator of equation (133) vanishes 
This leads immediately to the requirement that equation ( 1 50) can be satisfied which is contrary to hypothe- 
sis. It remains only to be shown that Q 2 (4> i ) is everywhere positive if there is no intersection. 

Again using the generic function f, let 


f = 7 cos <J>j + q - cos(<J>i + 0) 


so that extrema of f occur when t = a, i.e.. 
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(7 - cos {)) 

COS <J>! = ± 


sin cj) 


sin d 
r 


Then the extrema of f are given by 


fox = ± r + q 


so that, from (133), the extrema of are 


Q2.CX — Qi 


e 1 (r q)(r + q) 
2-y( ± r + q) . 


Suppose, first, that 

Q,„ = Q, [. + ^ ] 


and if q = r + e 2 , then 


Q2,cx 




— Qi 


as e becomes small, while if — q — r + rr then 


Q2. 


C\ 


= Q. 




as s becomes small. 


e,(Q,-Q.Q/(Q3e3)1 

(e,Qi)/(Q 3 e 3 ) J ^ 
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Next, let 


Q:.cx — Qi 


e i(r + q) 
2 7 


and q = r + e 2 then 


Q2 


Q, 


e 1 ( 2q — e 2 ) 


2y 


Q, 


[-•f] -»[' 


e » (Q 1 QjjAQjg 
(Qiei)/(Q 3 e 3 ) 


■*]- 


Q.i 


as e becomes small, while if — q = r + s 2 then 


Q2.ex 



Qi 


as £ becomes small. 

The conclusion is, then, that if the orbits do not intersect, Q 2 is bounded between the (positive) 
values of Q, and Q,. Since non-intersecting orbits have no singularities in the e 2 (cf>i) function and positive 
Q 2 ((t>,) values, we conclude that tangential transfer is possible between any two non-intersecting orbits. 
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